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Abstract. The isotropization process of a collisionless plasma with an electron temperature anisotropy along 
an external magnetic field (r| S> Tl, \\ and 1. with respect to the background magnetic field) and isotropic 
protons is investigated using a particle-in-cell (PIC) code. Restricting wave growth mainly parallel to the external 
magnetic field, the isotropization mechanism is identified to be the Electron Firehose Instability (EFI). The free 
energy in the electrons is first transformed into left-hand circularly polarized transverse low-frequency waves 
by a non-resonant interaction. Fast electrons can then be scattered towards higher perpendicular velocities by 
gyroresonance, leading finally to a complete isotropization of the velocity distribution. During this phase of the 
instability, Langmuir waves are generated which may lead to the emission of radio waves. A large fraction of the 
protons is resonant with the left-hand polarized electromagnetic waves, creating a proton temperature anisotropy 
Til' < Tj^. The parameters of the simulated plasma are chosen compatible to solar fiare conditions. The results 
indicate the significance of this mechanism in the particle acceleration context: The EFI limits the anisotropy of 
the electron velocity distribution, and thus provides the necessary condition for further acceleration. It enhances 
the pitch-angle of the electrons and heats the ions. 
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1. Introduction 

The energization of large numbers of particles to high en- 
ergies is still an unsolved problem in solar flare research. 
The vast diversity of observational data imposes severe 
constraints, for which any acceleration model has to ac- 
count. The current models can broadly be assigned to 
three different categories: Shock acceleration, acceleration 
by parallel electrical fields and stochastic acceleration by 
cascading MHD waves. For a review of these mechanisms, 



see e.g. [Miller ct al (1997] ). 

One restriction on the acceleration model is imposed 
by the time scales at which the process has to operate: 
Hard X-ray observations revealed spiky structures as short 
as « 400 ms ( |Kiplinger et al, 1984| ; pVIachado et al, 1993| ). 
Assuming a thick-target emission process for the radia- 
tion, the electrons have to be accelerated to energies above 
20 keV within the given time in each of these spikes. It is 
believed that these 'energy release fragments' are the ba- 



sic constituent of the entire impulsive phase ( Machado el 



al, 1993). Further fragmentation of the energy release and 
electron acceleration process was proposed based on deci- 



metric radio emission (Benz, 1985), with temporal struc- 
tures well below « 10 ms. 
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Assuming an isotropic electron velocity distribution. 
Miller et al (1996) showed that stochastic acceleration can 
account for energization of electrons from thermal ener- 
gies up and above 20keV in less than 400 ms. However, 
the actual particle acceleration mechanism in this model, 
Transit-Time Damping (Fisk, 1976; 3tix, 1992), increases 
only the velocity component parallel to the external mag- 
netic field. In absence of an additional pitch-angle scatter- 
ing mechanism, this would lead to an anisotropic electron 
velocity distribution, reducing the efficiency of the accel- 
eration mechanism. Due to the low density and the high 
temperature, the mean free path between Coulomb colli- 
sions in a solar flare plasma is large. The scattering mech- 
anism has therefore to be of a different nature, presumably 
a wave-particle interaction. 

Predominant acceleration along an external magnetic 
field Bq compared to the acceleration in perpendicular 
direction is not only a problem of Transit-Time damp- 
ing, but common to all the aforementioned acceleration 
mechanisms. The isotropization in case of particle accel- 
eration by a sub-Dreicer DC field has be en investigated by 
Moghaddam-Taahcri and Gocrtz (1990 ). Assuming quasi- 
linear diff usion and a wave s pectrum due to the runaway 
electrons ( Candy et al, 1983 ), they found that the actual 
isotropization mechanism is anomalous Doppler resonance 
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of fast particles and that a bump in the reduced velocity 
distribution is formed. 

In this paper, the isotropization process of an elec- 
tron temperature anisotropic plasma without further en- 
ergy input is addressed. The free energy in the anisotropic 
electron velocity distribution may drive an instability. The 
instability would lead to wave growth which in turn can 
scatter the particles in velocity space. One proposed in- 
stability for a electron temperature anisotropic plasma 
is the Electron Firehose Instabilit y (EFI) ([Hollweg and 



Volk, l [970| ; [Paesold and Benz, 1999| ; |Li and Habbal, 2000| ) 
It is a kinetic version of the MHD firehose instability 



( Rosenbluth, 1956 ; Parker, 195S ) with a lower threshold 
in electron velocity anisotropy. But is this mechanism ac- 
tually responsible for the isotropization of the electron 
velocity anisotropy? The previous investigations of this 
mechanism were based on linear theory, thus stopped be- 
fore wave growth effects the particle distribution. If the 
EFI actually grows or if another instability isotropizes the 
plasma quicker cannot be answered within the framework 
of linear theory. Also questions about wave growth satu- 
ration are beyond reach of linear theory. 

Here, the isotropization mechanism of an electron tem- 
perature anisotropy is investigated by means of particle- 
in-cell (PIC) simulations, giving access to the non-linear 
effects on the particle distribution as well as the satura- 
tion of the instability. The parameters for the simulation 
are chosen to be compatible with parameters expected un- 
der solar flare conditions. While rough estimates for den- 
sity, ambient magnetic field and electron temperature ex- 



ist (e.g. Miller et al, 1997), values of the actual anisotropy 
depend on the exact details of the acceleration mechanism. 
The anisotropy chosen here is T£/Tj^ = 20 in accordance 



with earlier estimates ( Paesold and Benz, 1999 ). The sim- 
ulated plasma corresponds to a population of electrons 
accelerated from 3keV to 25keV without any instabilities. 

The paper is organized as follows: Section || de- 
scribes the simulation model, followed by a presentation 
of the overall development of the isotropization process 
in Sect. ^. The responsible instability is then identified 
by comparison with linear theory in Sect. ^. Non- linear 
processes are discussed in Sect. ||, followed by an investi- 
gation of the energy flow in Sect. ^. A discussion of the 
applicability of result onto real plasmas (Sect. ^ and of 
the relevance of the described mechanism summarizes the 
results in Sect. 0. 



2. Simulation l\/lodel 



The simulations were performed with the fully 3D 
relativistic electromagnetic particle- in-cell code par-T 
(iMessmer, 2000|) a parallel implementation of TRISTAN 
( Buncman, 1993 ). As both particle species are expected to 
play a role in the electron firehose instability, they both 
have to be treated kinetically. This prevents the use of 
hybrid simulations which treat one particle species as a 
fluid. 



A PIC code traces the trajectories of a represen- 
tative number of particles under the influence of the 
Newton-Lorentz force in their self-consistent electromag- 
netic fields. The particles can be placed anywhere within 
the computational domain, whereas the field quantities are 
located on a regular spatial grid. Subgrid resolution, e.g. 
in order to determine the forces at the particle positions, 
is obtained by linear interpolation. For an introduction on 



PIC see e.g. Birdsall and Langdon (1985) 



Time integration for both particles and fields is per- 
formed by a leap-frog scheme. The fields are staggered 
in space { Yee, 1966 ), allowing e.g. to update the B-field 
by computing the curl of the E-field, using finite differ- 
ences. This method has the advantage of being simple to 
implement and additionally to keep VB = 0, if it was sat- 
isfied initially. In oder to update the E-field, the current 
density has to be determined in addition to the curl of 
B. Instead of solving the Poission equation explicitly, the 
herein used code applies a charge conserving current de- 
position algorithm ( Villasenor and Buneman, 1992| ). This 
scheme allows to update the E-ficld entirely from informa- 
tion within a few grid cells. Additionally it keeps V E — p 
satisfied, if it was satisfied initially, down to machine ac- 
curacy. 

The size of the simulation box is Lj. x Ly x Lz = 
8 A X 8 A X 768 A, where Lj is the length of the system 
in dimension j, A = 0.13c/we represents the cell size, 
oje is the electron plasma frequency and c is the speed 
of light. Periodic boundary conditions in all three dimen- 
sions are applied on both particles and fields. The total 
particle number is 2 • rip = 3.17 • 10® particles, an aver- 
age of 32 electron-proton pairs per cell. The particles are 
placed uniformly within the computational domain. The 
simulation time step is WgAi — 27r/100. Along the z-axis, 
the external magnetic field Bq is applied, leading to a ra- 
tio flg/uje = 0.13, where fie is the electron cyclotron fre- 
quency. Assuming a particle density of n — 5 ■ lO^'^cm"^, 
the external magnetic field is Bo — 91 Gauss. As the 
isotropization is expected to take place on proton time 
scales, an artificial proton/electron-mass ratio mp/nie = 
49 is chosen, a tribute to computing time. The proton 
plasma frequency is therefore ojp/uje = 0.14 and the pro- 
ton cyclotron frequency flp/cue = 0.0026. The initial elec- 
tric field is E'o = everywhere within the computational 
domain by placing electrons and protons at the same po- 
sitions. 

Initially, the electron temperature parallel to the ex- 
ternal field is T| = 2.3-10^ K while the perpendicular elec- 
tron temperature T|[ is identical to the isotropic proton 
temperature — 1.2 • 10^ K. The electron temperature 
ratio is therefore T^/T^ w 20. The corresponding thermal 
velocities are for electrons vpc = 0.2 and w^/c = 0.045. 
The Debye lengths are A^iji = 1.53A and Xd,± = 0.34A. 

By choosing a rod shaped geometry, wave numbers 
at oblique angles are limited to values k > /c™'" — 
k^™Lz/ = 42wp/c, where fc™'" is the minimum wave 
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number in parallel direction. The simulation can therefore 
be considered as a ID simulation. 

The instability is expected to grow on timescales of 
proton cyclotron periods. This makes it necessary to run 
the simulation at least for several proton cyclotron peri- 
ods. On the other hand, the free energy is carried by the 
electrons, requiring resolution of the electron plasma fre- 
quency. Even with an artificially low mass ratio, this leads 
to several thousand timesteps per inverse proton cyclotron 
frequency. To overcome the enormous computing effort, a 
parallelized simulation code was applied. Ideally, such a 
code requires a computing time Tp — r/Np, where Np is 
the number of processors, as compared to r, the time re- 
quired by the sequential code. In a more realistic model, 
the speedup S = r/Tp < Np, mainly due to the overhead 
introduced by the communication among the processors. 
However, in case of a par ticle-in-cell code with charge con- 
serving current update ( Villasenor and Buneman, 1992 ) 
only local information is used to solve the field equations. 
Consequently most of the time is spent in pushing particles 
and updating field quantities, letting communication be of 
minor influence. The code does therefore not benefit from 
a highly sophisticated interconnection network. It can eas- 
ily be run on a cluster of workstations, which are usually 
much cheaper and more available than commercial paral- 
lel computers. The simulations were therefore performed 
on a Beowulf-Cluster featuring 48 Pentium-Ill processors. 
The total simulation time of 120'000 time steps requires 
24 hours, including diagnostics. At a lower limit of the 
parallel efficiency of e « 0.8 on 48 processors ( Messmer 
2000|)7|this corresponds to about 920 hours of the original 
TRISTAN code on a single Pentium III processor. 



3. Overall Development 

Figure |l| shows the development of the local temperatures 
for both the electrons and the protons during the instabil- 
ity. As both the electron and the proton velocity distribu- 
tion may not stay Maxwellian throughout the instability, 
a global temperature is not defined. The local tempera- 
ture is defined as the standard deviation in one velocity 
component for a set of particles in a small subdomain of 
the computational domain. The temperature of the distri- 
bution is then defined as the average of all local tempera- 
tures. 

During the simulated time of 18^1^^, the inital elec- 
tron temperature ratio reduces to T^/Tj^ < 1.2. The 
free energy in the electrons goes partially into the waves, 
but mainly into the perpendicular electron temperature. 
This will be shown in Sect. ^. Additionally the protons 
are slightly heated in all directions. At about 80~^, the 
proton velocity distribution is becoming anisotropic with 

> T^. The oscillations are an indicator of energy trans- 
fer between the protons and the waves. 

For comparison, the parallel temperature development 
of an isotropic plasma is also plotted in Fig. |l|. As 
expected, the temperature does practically not change 
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Fig. 1. Temporal development of the local electron tem- 
peratures (top) and proton temperatures (bottom), com- 
pared to the temperature development in case of an 
isotropic plasma (thin dash-dotted lines). Temperatures 
are shown along the three spatial dimensions x, y, z (dot- 
ted, dashed, sohd). 



throughout the simulated time. This is a first indicator 
of energy conservation within the simulation. 

Figure ^ shows the time history of the perpendicular 
magnetic field energy Bj^ = JJJ{Bl + By)dx dy dz. The 
development of B^ for a simulated plasma with isotropic 
temperatures r| = = but otherwise identical pa- 
rameters is again included for comparison. 

After a slow increase during 857"-'^, the energy in B± 
increases explosively, saturates at 12. 5$!"^ and then starts 
to oscillate, indicating the growth of transverse (elec- 
tro)magnetic waves. 



4. Comparison with Linear Theory 

In order to investigate the isotropization further, the ex- 
cited waves have to be identified. The wave numbers are 
determined by Fourier transform of B^, one component 
of the perpendicular B-field. Figure ^ shows the temporal 
evolution of the spatial mode kc/ujp — 1.32 propagating 
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Fig. 2. Temporal development of -B^/Bjj for the 
anisotropic electron velocity distribution with T||^/T^ w 20 
(solid) and = T^, compared to the development of 
for an isotropic electron distribution with T^^ = 
Tl = TP (dotted). 
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Fig. 3. Temporal development of \Bx\ for the wave mmr- 
ber kc/uip — 1.32. 



parallel to Bq. After an initial phase of w Iflp ^, a regular 
oscillation starts with exponentially growing amplitude. 
Close to 13rJ-\ the wave growth stops and the oscillation 
continues. 

To identify the wave mode responsible for the growing 
field J5_L, the growth rates of different wave numbers are 



compared to the values obtained by WHAMP ( Ronnmark. 
1982()7]rhis program allows to evaluate the plasma disper- 
sion function numerically for arbitrary temperature ratios. 

Figure ^ shows a comparison of the growth rates be- 
tween linear theory and PIC simulation. The linear theory 
results are determined for the initial plasma parameters. 
The growth rates of the PIC simulation are measured by a 
linear fit in the time series of spatial Fourier transforms of 
Bx in the time interval 2 — Sfi"^. The three growing waves 



Fig. 4. Growth rate of B^ determined from simulations 
with mp/nie — 49 (asterisk) and mp/rrif, = 81 (boxes) 
compared to results obtained from linear theory for 
mp/nie = [49,81,1836] (dotted, dashed, dash-dotted) for 
the initial electron temperature ratio Tn/T'[ « 20. In ad- 



dition, the growth rate for T|/TJ 
is shown (dotted, thin). 



6 and mp/m^ 



49 



and the suppressed growth at kc/ojp = 1.76 are well in ac- 
cordance with linear theory. However, the growth rates 
for kc/ojp = [0.44,0.88] are too low, whereas the one for 
kc/uOp = 1.32 is too high. This will be investigated in 
Sect. 



5.1 



In order to check the influence of the artificial mass 
ratio, the growth rates are determined for an increased 
mass ratio of nip/me = 81 and otherwise identical plasma 
parameters (see Fig. ||). Due to the larger mass ratio and 
the fixed simulation box size of lOOc/we, the resolution in 
wave numbers decreases. Again, the waves grow in accor- 
dance with linear theory, but the determined growth rates 
are lower than expected for kc/ujp = 0.59 and a bit larger 
for kc/ujp = 1.18. For comparison, the growth rates for the 
real mass ratio mp/m^ = 1836 is included in Fig. ^. The 
hypothetical light weight protons inhibit wave growth at 
large k. 

Figure ^ shows a comparison of the dispersion rela- 
tion determined both by WHAMP and from the simula- 
tion. Due to the limited resolution in k-space, only 3 sam- 
ples lie within the range of non-zero growth (see Fig. 
According to linear theory, the period of the sample with 
smallest wave number, kcjujp = 0.44, is 4517"^ and thus 
much longer than the saturation time of the instability. 
The frequency can therefore not be measured in the sim- 
ulation and only an upper bound can be given. For the 
other two samples, the frequencies are determined by re- 
moving the exponential growth and measuring the time 
between the first two amplitude minima, thus the first 
half period. Like the growth rates, the wave frequencies 
agree well with the values determined from linear theory. 
However, they are both too low. 
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Fig. 5. Comparison of the dispersion relations obtained 
by the simulation (asterisk) and by WHAMP (dashed). 
For comparison, the dispersion relation for the real mass 
ratio mp/nif, — 1836 is included (dash-dotted). 

Table 1. Comparison of growth rates [70.44,70.88,71.32] 
between simulation (bold) and linear theory at kc/ujp — 
[0.44, 0.88, 1.32] and wave number of maximum growth km 
for different electron temperature ratios T^/T^ . The other 
plasma parameters are identical to Sect. 
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Additional simulations with varying initial tempera- 
ture ratio have been carried out to test the agreement with 
linear theory. Table |l| shows the growth rates for varying 
temperature ratios T^^/T^, but otherwise identical param- 
eters as in Sect. ^. A larger anisotropy leads to a shift of 
the maximum growth rate towards smaller wave numbers 
(see also Fig. ^). 

The good agreement of both dispersion relation and 
growth rate indicates that the instability under investiga- 
tion is indeed the Electron Firehose Instability. 



5. Non-Linear Effects 

5.1. Decreasing Temperature Anisotropy 

The previous section showed that the rough estimates for 
the growth rates and the wave frequencies agree well with 
linear theory. The discrepancies can be qualitatively un- 
derstood as an effect of the temperature anisotropy reduc- 
tion: As mentioned in the last section, the growth rates 
are determined by linear fit between 2 — 8f2~^. They are 
therefore average growth rates over the whole time in- 
terval. However, during this time the temperature ratio 
is reduced to a value T^/T^ « 6 (see Fig. This has 



Fig. 6. Temporal development of the wave frequency for 
kc/uip — 1.32 compared to the frequencies expected from 
linear theory according to the instantaneous temperature 
anisotropy. 

an influence on the growth rates of the EFI, shown in 
Fig. IJ: For lower anisotropy, the maximum growth rate 
is shifted towards larger wave numbers. For the reduced 
anisotropy, the growth rates for kc/cup < 1.1 are lower and 
for kc/ujp > 1.1 are larger than for the initial anisotropy. 
Averaging the growth rates at kc/uop < 1.1 thus underes- 
timates the initial growth rates, whereas for larger wave 
numbers the growth rates are overestimated. This agrees 
with the measurements in Fig. ^ and Tab. |l|. 

The reduction of the anisotropy has not only an ef- 
fect on the growth rate, but also on the real part of the 
wave frequency. Fig. ^ includes the dispersion relation for 
the reduced temperature ratio Tu/T'[ w 6, showing a fre- 
quency decrease for all wave numbers. Measuring the time 
between two minima in spectral power will only yield an 
average frequency. 

Figure ^ shows the development of the average fre- 
quency for the wave number kc/ujp = 1.32 (see also Fig. 
The horizontal bars indicate the time interval between two 
amplitude minima. The errorbars in frequency represent 
the uncertainty in determing the time of the minimum 
amplitude. The average frequencies agree well with the 
frequencies expected from the linear theory for the instan- 
taneous temperature anisotropics. 

5.2. Saturation 



Based on linear theory, HoUweg and Volk (197G| ) give 
an approximate instability criterion for growth of a left 
hand polarized electro-magnetic wave in a temperature 
anisotropic plasma. 



< 



(1) 



where /3|j — 2(ijJe/^^e)^ • {vf\/c)'^ is the ratio of thermal 
to magnetic pressure, and = I — T^/T^ is the elec- 
tron temperature anisotropy. This criterion is strictly valid 
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Fig. 7. Temporal development of 1 — P^^Ae for nip/me — 
49 (solid) and mp/nie = 81 (dashed). The dotted hne 
represents the instability criterion given by Eq. M). 
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Fig. 9. Number of particles in a velocity range of 
Aw||/c = 0.04 centered at a fixed W||/c = [0, 0.16, 0.30, 0.50] 
[solid, dotted, dashed, dash-dotted] as a function of time. 



only for small anisotropics |1 — /3|jAe| ^ 1. What happens 
in case of a larger anisotropy? The threshold is identi- 
cal to the MHD firehose instability criterion in case of an 
isotropic proton velocity distribution. The simple physi- 
cal interpretation is that the centrifugal force exerted by 
the electrons moving along a disturbed magnetic field line, 
destabilizes the wave. If the anisotropy is not large enough, 
the magnetic field stress and the thermal pressure perpen- 
dicular to the field line act as restoring forces. If the ve- 
locity anisotropy is large and the magnetic field is small 
enough, the restoring force is too weak and thus the wave 
destabilizes. 

Figure ^ shows the temporal development of the 
anisotropy and instability criterion given in Eq. (|^), as- 
suming a constant background magnetic field. The initial 
setup lies within the unstable region, triggering the in- 
stability. As soon as the stable region is reached at time 
12Q~^, wave growth stops. The instability criterion then 
states, that the anisotropy in the electron velocity distri- 
bution has mostly been eroded. 



5.3. Electron Velocity Distribution 

An advantage of self-consistent kinetic simulations is the 
possibility to investigate effects on the particle distribu- 
tions due to wave growth. 

Figure || shows the electron velocity distributions at 
times t = [0,6.5,7.8,10.7, 12.7, 18.0]O-^ It demonstrates 
that the initial velocity distribution isotropizes in the 
course of the simulated 18fl~^. This is a result of the 
EFI growing at the same time scales as the isotropiza- 
tion takes place. A closer look at the velocity distribution 
reveals two distinct features: on one hand, for the bulk of 
the electron velocity distribution, v± increases at small W|| . 
On the other hand, particles with |?;||| 3> are ejected 
from the distribution to larger pitch angles in narrow ve- 
locity ranges (see Fig. ^, tflp — 7.8). The first feature is 



a sign of non-resonant wave-particle interaction, whereas 
the second is a sign of resonant wave-particle interaction. 

5.3.1. Non-Resonant Wave-Particle Interaction 

According to HoUweg and Volk (1970), most of the elec- 
trons gain perpendicular momentum by a non-resonant 
wave-particle interaction. The effect of such an interaction 
is pitch-angle scattering of electrons with no preferred v\\ , 
increasing at the expense of |v||| by randomly fluctu- 
ating wave fields. The non-resonant pitch-angle scattering 
can be seen in Fig. ^before t « 7.8r2~^, where the perpen- 
dicular velocity increases at all ?;|| concurrently. But is it 
actually pitch-angle scattering or heating in perpendicular 
direction? 

Figure |^ shows the temporal development of the elec- 
tron number density iV in a velocity range of Au|| /c — 0.04 
centered at a fixed v\\/c — [0,0.16,0.30,0.50]. The par- 
ticle numbers are normalized to the particle numbers at 
t — 0. The particle number for W||/c = 0.16 « ^f]/'' remains 
about constant throughout the growth of the instability. 
The particle number at low velocities (w||/c = 0) increases 
with time, while the particle numbers for W||/c 3> v^^ de- 
crease. At W||/c — 0.50 an increase of particles is visible 
at 8ft~^ due to the particles resonant with the wave. This 



will be discussed in Sect. 5.3.2| . 

The increase of particle number at low |t;||| and the 
decrease at high |u||| indicates that the particles are not 
simply heated in perpendicular direction, but that they 
are scattered towards lower parallel velocities. How is the 
magnetic field built up? The scattering of a single particle 
due to randomly fluctuating fields corresponds to an accel- 
eration, inducing a magnetic field pulse, which propagates 
as a wave. Most of these waves are not linear modes of the 
given plasma and are therefore highly damped. However, 
some of these waves correspond to the linear eigenmodes 
of the plasma with positive growth rates. After a certain 
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Fig. 8. Projection of the electron velocity distribution along Vy at times t = [0, 6.5, 7.8, 10.7, 12.7, 18.0]J1~^ (left to 
right, top to bottom). The contour levels are at [1, 10, 100, 1000] particles. The ambient magnetic field points from left 
to right. Overplotted are the resonance curves for the three modes kc/up = [0.89, 1.32, 1.76] (dotted, solid, dashed). 
The aspect ratio is chosen such that circular contour lines correspond to a thermal distribution. 



build-up time, the amplitude of the induced wave is large 
enough to influence the particle velocity distribution. In 

the present simulation, this happens at about 7.5f7p ^. At 
this time, fast electrons start to gyroresonate with the 
growing wave. 



5.3.2. Resonant Wave-Particle Interaction 

In presence of an external magnetic field, anomalous gy- 
roresonance of the electrons in the form 

w — kv\\ = —s— (2) 
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can be expected, where uj, and k are frequency and wave 
number of the resonant wave, W|| is the parallel velocity of 
the resonant electrons, s is the harmonic resonance num- 
ber, i7e the non-relativistic electron cyclotron frequency 
due to the external magnetic field and 7 the Lorentz fac- 
tor of the particle. 

The frequencies for the fastest growing modes are de- 
termined by linear theory (see Sect. The relativistic 
gyroresonance condition can be solved numerically for an 
electron with parallel speed /3 = W||/c, 

$ = Kf3^ + 1 (3) 
7 

where $ ~ uj/V,e and k ~ kc/fle ( ^teinacker and Miller, 
1992|riFor a detailed discussion of gyroresonance see e.g. 
Benz (1993[) . 

The resulting resonance curves for the growing waves 
due to the EFL kc/ujp = [0.88,1.32,1.76] are plotted in 
Fig. ||. Although the wave at kc/uip = 1.76 does not grow 
initially, it becomes unstable as soon as the temperature 
ratio is small enough. This is the case at about 8.3rip^. 
The resonance curves are also affected by the reduction in 
temperature anisotropy, but the effect is small. 

Up to « 6.5ilp^, no specific features in the velocity 
distribution can be seen close to the resonance curves. 
B± is growing, but it is not yet strong enough to affect 
the electron velocity distribution. At « 7.851"^, particles 
with velocities close to the resonance velocity of the fastest 
growing wave kc/cup = 1.32 are in resonance and as a con- 
sequence, are scattered towards higher perpendicular ve- 
locities. According to linear theory, the growing wave is a 
left hand polarized electromagnetic wave. The interaction 
with the electrons is therefore anomalous Doppler reso- 
nance. At a later stage, the wave amplitude of the slower 
growing waves can also interact with the electrons, scat- 
tering them to even lower and higher vj_. 

Only the high velocity tails of the distribution are af- 
fected by the anomalous Doppler resonance. The major 
fraction of the electron particle distribution gains perpen- 
dicular velocity by the non-resonant interaction. 

During the phase of anomalous Doppler resonance, 
the reduced electron velocity distribution can generate a 
bump on the velocity distribution tail. Figure |l^ shows 
the temporal development of the reduced electron veloc- 
ity distribution parallel to the external magnetic field, 
fBim^t) = / f{v±,vii,t)dv±, where f{v±,vii,t) is the 
electron velocity distribution at time t, between times 
S.im-'^ and 9.08^-1 at intervals of 0.1617"^ for high ve- 
locities. Throughout this time interval, a plateau persists 
in the distribution between the two resonance velocities for 
kc/ujp = [0.88,1.32]. At the time SASfl'^, a clear bump 
between wy/c = —0.62 and fy/c = —0.66 is visible in the 
reduced electron velocity distribution. Due to the symme- 
try of the problem, a similar plateau with bump exists for 
the positive velocities 0.62 < W||/c < 0.66. 

A bump on the reduced velocity distribution leads to 
the bump-on-tail instability, generating Langmuir waves. 
The wave numbers, fc^, of the Langmuir waves generated 




Fig. 10. Evolution of the reduced electron velocity dis- 
tribution /_r(^^||) between t — 8.1157"^ (bottom) and 
t = 9.03fJ-^ (top) at intervals of 0.16^-^. The bold line 
indicates the distribution at t = 8.4317^ ^ with a clearly 
visible bump. The vertical lines represent the resonant ve- 
locities v\\ and < v±/c < 0.4 for the wave numbers 
kc/ujp = 0.88 (dashed) and kc/ujp = 1.32 (dotted). 

by particles with velocity wy can be estimated by combin- 
ing the Langmuir dispersion relation with the Cerenkov 
resonance condition. 



Figure |Il| shows the temporal evolution of the paral- 
lel electric field amplitude Ez for different wave numbers. 
Only field amplitudes exceeding the mean amplitude by 4cr 
are shown, a being the standard deviation. Overplotted 
are the Langmuir wave numbers for the beam velocities 
W||/c = 0.62 and v\\/c = 0.66. The bump-on-tail distribu- 
tion in Fig. |l^ leads obviously to the emission of Langmuir 
waves. 

One could expect these longitudinal waves (L) to cou- 
ple to the transverse EFI waves (T), generating observ- 
able transverse radio emission, L + T ^ T. However, the 
large wave numbers of the Langmuir waves and the low 
frequencies of the EFI waves make it impossible to sat- 
isfy the parametric equations for observable radio waves. 
Whether these waves can couple to other plasma waves in 
a more complicated way, yielding narrowband radio emis- 
sion, needs to be investigated further. 

5.4. Proton Velocity Distribution 

The initial proton velocity distribution is assumed to be 
isotropic at a temperature of = T'[. Although protons 
do not contribute any free energy to the instability, they 
play a significant role in the course of the instability. As 
the growing wave due to the EFI is left hand polarized, the 
protons come in normal Doppler resonance with the wave. 
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Fig. 11. Temporal evolution of the parallel electric field 
amplitude for different wave numbers. Dark corre- 
sponds to locations where the amplitude exceeds the aver- 
age field amplitude by 4(7. Overplotted are the wave num- 
bers for Langmuir waves generated by a parallel propagat- 
ing weak beam with W||/c = 0.62 (dotted) and v\\/c = 0.66 
(dashed). 




Fig. 12. Projection of the protons in phase space to the 
plane Vx vs. z aX t ^ lO.Tfip^. Contour levels are at 
[1, 10, 100] particles. The ambient magnetic field points 
from left to right. 



Unlike electrons, where only the high velocity tails of the 
distribution are affected by resonance, the resonance ve- 
locity for the protons lies within the bulk of the proton ve- 
locity distribution. Therefore the protons are expected to 
absorb energy from the wave. The non-coUisional electron- 
proton coupling leads to a partial transfer of thermal en- 
ergy from the electrons to the protons. This energy tran s- 
fer was already conjectured by HoUweg and Volk (197C ). 



Due to the low thermal velocities of the protons, 
the gyroresonance velocity can be estimated in a non- 
relativistic limit, simplifying Eq. (j^) to 



UJ - 



kc/uip 



-CA 



(5) 



for resonance with an left hand polarized wave, where 
ca = c^lp/iUp is the Alfven-speed. The initially 
fastest growing waves {kc/ujp;uj/V,p) = (0.88; 0.68) and 
(1.32; 1.82) have resonance velocities v^^/v^ ~ —1.0 and 
1.8. Both waves could be in resonance with some of the 
protons but due to the small amplitude of the waves, their 
influence is small. 

Due to the isotropization of the electron velocity distri- 
bution, the wave frequencies at fixed wave number change 
and therefore the proton resonance velocities have changed 
to vii/vP = -2.1 for (0.88; 0.33) and to wy/vP = 0.2 for 
(1.32; 1.07) at time 7.3^l~^. The resonance velocity for 
kc/uip = 0.88 is already far out in the proton velocity 
distribution and does not affect the proton velocity distri- 
bution any more. For the wave (1.32; 1.07) on the other 
hand, a large fraction of the proton distribution can be in 
resonance. 



Figure 12 shows the proton distribution in a projection 
of phase space at 10.7Qp^. The spatial wave pattern cor- 
responds to kc/up = 0.88, which is not in resonance with 
the protons. The wave with kc/ujp = 1.32 on the other 
hand is in resonance with the bulk of the proton veloc- 
ity distribution, leading to pitch angle scattering of the 
protons towards larger \v±_\ while reducing 

6. Energy and Temperature Development 

The transfer of free energy from the electrons to the pro- 
tons through wave-particle interaction can be observed by 
considering the relevant energies in the simulated system. 

Figure |l3| shows the development of the different ener- 
gies per particle, demonstrating the transfer of free elec- 
tron energy into waves and from there into proton kinetic 
energy. The kinetic energy of species i is given by 



El ^ —Y^mpc^i^, 



1), 



(6) 



where 7-/ is the Lorentz- factor of particle j. The electro- 
magnetic field energy is given by 



Ef 



1 

2n^ 



[E"^ + c^B^)dxdy dz, 



(7) 



where E^ = E^ + E^ + E^ and 3^=31+3^ + B^. 

First it should be noted that throughout the simulated 
time the total energy is conserved to less than 2%. 

Initially, an electron has on average a total energy of 
« 25 kcV, with 23 keV in parallel direction and 1 keV 
in each perpendicular direction. The protons are initially 
isotropic with 1 keV per proton in all directions. 
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Fig. 13. Energy development per particle: Total energy 
Ef + E^ + E^ (dotted), electron kinetic energy (solid), 
electromagnetic field energy E^ (dash-dotted), proton ki- 
netic energy E^ (dashed), and E^ + E^ (dashed, thin). 



Up to « lOrip^, the electron kinetic energy remains 
about constant. At the same time the parallel electron 
temperature is reduced significantly, thus most of the free 
energy is used to heat the electrons in perpendicular di- 
rection. However, some of the electron energy is converted 
into magnetic field energy. Due to the high mobility of 
the electrons, large scale electric fields do not build up. 
Between 10 — 12n~^, energy goes mainly into the mag- 
netic field. After saturation, the electron energy remains 
constant, indicating that the electrons play a passive role 
in the later development. This is also supported by the 
almost constant E^ + E^ . 

The protons do not increase their energy significantly 
up to about lOilp ^. At that time, the bulk of the protons 
becomes resonant with the waves, increasing their kinetic 
energy at the cost of the wave. However, due to the reso- 
nant character of the interaction between protons and the 
waves, the proton kinetic energy can also be fed back into 
the magnetic field, leading to the oscillations seen after 



15. 



7. Real Plasmas 

The isotropization process of an electron temperature 
anisotropy has been investigated in case of a plasma 
with mp/rrie — 49. What happens in a real plasma with 
mp/rrie — 1836? As shown in Fig. ^ and Fig. |[ the ef- 
fect of an increased proton mass is wave growth at larger 
k. Additionally the maximum growth rate will be slightly 
larger. 

The shift of the maximum growing wave number for re- 
duced anisotropy, as well as the reduction of the wave fre- 
quencies at fixed wave numbers, will also occur in a plasma 
with real mass ratio. E.g. a reduction of the electron tem- 
perature ratio from T^/Tl « 20 to T,f/Tl « 6 will shift 



the wave number of maximum growth from kc/ujp = 1.2 
to kc/ujp — 1.8. 

What happens to the protons? Keeping the ratio 
U,e/uje constant, the Alfven-speed ca has to be scaled by 
^mj., where is the ratio between artificial and real pro- 
ton mass. At the same time, the thermal velocity has to 
be scaled by ^Jm^, thus the ratio v\\/v^ remains constant. 
However, for a fixed wave number, the frequency decreases 
due to the different dispersion relation. It was shown in 



Sect. 5.4, that the fastest growing wave was in resonance 



with protons much slower than the thermal velocity. This 
is still true in case of a plasma with real mass ratio, e.g. 
for a temperature ratio of T^/T^ k, 6, the fastest grow- 
ing mode is {kc/uip\Lo /Vlp) ~ (1.7; 0.9) and the resonance 
velocity is v^/vP — 0.16. 

Due to the resonant absorption of the waves by the 
protons and the broader wave spectrum excited in a real 
plasma, the heating of the protons due to the EFI may 
be larger than predicted by the simulations. Additionally 
it was demonstrated in the previous section, that the ab- 
sorption of field energy by the protons has only minor 
influence on the electrons. As the only drain of free elec- 
tron energy is either isotropization or wave growth, and a 
broader spectrum of waves is excited in a real plasma, the 
simulated electron temperature isotropization time can be 
assumed to be an upper limit. Faster isotropization is 
confirmed by a simulation with increased mass ratio of 
nip I = 81 (see Fig. 0). 

8. Summary 

The preferred acceleration along an external magnetic 
field compared to the perpendicular direction causes an 
unsolved problem in solar flare particle acceleration mech- 
anisms. Without additional scattering mechanisms, this 
leads to an anisotropic electron velocity distribution, 
which can in turn reduce the efficiency of the accelerator. 
The isotropization time of an anisotropic velocity distri- 
bution has to be short compared to the total acceleration 
time which is 0.4 s or less to energize the electrons to 20 
keV. 

Based on linear theory, the Electron Firehose 
Instability (EFI) has been proposed as a candidate for ve- 
locity space particle scattering. In the previous sections, 
the temporal development of the EFI and its influence on 
the particle velocity distribution has been shown, based 
on self-consistent electromagnetic PIC simulations. The 
plasma parameters were chosen to be comparable to ex- 
pected values in a solar flare plasma. 

In case of a rod shaped geometry, which allows long- 
wavelength modes to grow mainly parallel to the external 
magnetic field, the process responsible for isotropization 
is identified to be the Electron Firehose Instability. The 
identification is made by comparison of the growth rate 
and the dispersion relation between simulations and linear 
theory. 

The EFI is driven by the bulk of the electron velocity 
distribution, and is thus non-resonant. Most of the parallel 
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electron energy is thereby converted into perpendicular 
electron energy. However, some free energy is converted 
into magnetic field energy. 

The growing waves are left hand polarized electromag- 
netic waves with frequencies close to Qp. As soon as the 
wave amplitudes are large enough, they start to scatter 
the electrons in velocity space by anomalous Doppler res- 
onance. In this phase, the electrons emit most of their 
free energy as waves by reducing quickly their anisotropy. 
Wave growth stops after the anisotropy has been reduced 
below the instability criterion. 

The strong waves generated by the instability mainly 
show up in the perpendicular magnetic field. Due to their 
left hand polarization, they are normal Doppler resonant 
with the protons. The resonance velocities of the waves 
generated by the EFI arc in the bulk of the proton veloc- 
ity distribution, and can thus be easily absorbed by the 
protons. This leads to perpendicular heating of the pro- 
tons. 

At the time of the anomalous Doppler resonance of the 
electrons with the wave, the electrons build up bumps in 
the reduced velocity distribution. They arc the source of 
Langmuir waves, which in turn may yield observable radio 
emission. 

Under solar flare conditions, the simulated plasma cor- 
responds to a plasma with electrons of 2 keV in perpen- 
dicular direction and 23keV parallel to the external field. 
According to the simulations, this distribution isotropizes 
within about 20Q~^, corresponding to 2-I0~^s. The result- 
ing isotropic energy of the electrons is about 7keV. The 
isotropization is fast compared to the timcscalcs needed to 
accelerate the whole electron distribution to 25kcV, which 
is expected to be on many orders of magnitude larger time 
scales. The electron anisotropy originating from accelera- 
tion can therefore not be maintained. 

Due to the large growth rates of the EFI at oblique 
angles, occurrence of the oblique firehose instability in a 
real plasma is even more likely. However, it is not yet clear 
how much energy is carried by those modes and what their 
influence on the particle distribution is. Additionally, due 
to the artificial mass ratio, the isotropization time may be 
overestimated. The given isotropization time due to the 
EFI can therefore be considered as an upper limit for an 
anisotropy to exist. Increasing computing resources will 
allow to tackle the problem in a slap shaped or even fully 
3D geometry with larger mass ratios. This is the topic of 
current investigations. 

In these simulations, most of the free energy is used 
to heat the electrons in perpendicular direction. Little en- 
ergy goes into the protons, berating them in perpendicu- 
lar direction and cooling them in parallel direction. The 
isotropization process can therefore not be used to explain 
bulk heating of the protons. 

For the assumed plasma parameters, the Electron 
Firehose Instability limits the anisotropy. The particle 
distribution can therefore be assumed to be isotropic 
throughout the acceleration process. 
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